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We revisit the simplest model for dark matter. In this context the dark matter candidate 
is a real scalar field which interacts with the Standard Model particles through the Higgs 
portal. We discuss the relic density constraints as well as the predictions for direct and 
indirect detection. The final state radiation processes are investigated in order to understand 
the visibility of the gamma lines from dark matter annihilation. We hnd two regions where 
one could observe the gamma lines at gamma-ray telescopes. We point out that the region 
where the dark matter mass is between 92 and 300 GeV can be tested in the near future at 
direct and indirect detection experiments. 
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I. INTRODUCTION 

The possibility to describe the properties of the dark matter (DM) in the Universe with a par¬ 
ticle candidate is very appealing. This idea has motivated the theory community to propose a vast 
number of dark matter candidates and today we have many experiments searching for these candi¬ 
dates. The traditional way to look for dark matter is through direct detection where one expects 
to see the recoil energy from the scattering between the dark matter candidate and nucleons [1], or 
from the scattering between the electrons and the dark matter. One also could see exotic signatures 
at the Large Hadron Collider (LHC) associated with missing energy due to the production of dark 
matter. However, since one cannot probe the dark matter lifetime at colliders, this latter possibility 
is perhaps not the most appealing one. See Refs. [2-5] for reviews on dark matter candidates and 
corresponding experimental searches. 

The annihilation of the dark matter in the galaxy into gamma rays can provide a very striking 
signal which can be used to determine the dark matter mass and understand the dark matter 
distribution in the galaxy. One expects more photons in the center of the galaxy and the dark 
matter profile dictates how many photons one could expect in other regions of the galaxy for a 
given value of the annihilation cross section. Since the dark matter candidate does not have electric 
charge, the dark matter annihilation into monochromatic photons occurs at loop level, and it could 
be very difficult to observe these lines due to the continuous spectrum. See Ref. [6] for a recent 
review on dark matter annihilation into gamma rays. 

In the simplest dark matter model one has only a real scalar field [7], which is stable due to 
the existence of a discrete symmetry. This model has only two parameters (relevant for the DM 
phenomenology) and one can have clear predictions for direct and indirect detection experiments. 
Since this is the minimal theory for dark matter one should investigate all the predictions to 
understand how to test this model in the near future. This model has been investigated by many 
groups [8-25]. However, only recently it has been pointed out [24] that one can observe the gamma 
lines from dark matter annihilation in this context due to the fact that the final state radiation (FSR) 
processes are suppressed in some regions of the parameter space. This is by far not generic given a 
dark matter model, as photons from tree-level processes tend to dominate the spectrum. 

In this article we revisit the singlet dark matter model investigating all current constraints from 
relic density, invisible Higgs decays, direct and indirect detection. Here we complete the study 
presented in Ref. [24]. Our main aim in this article is to present the scenarios that can be tested 
using direct and indirect detection experiments. We focus on the discussion of gamma-ray lines and 
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point out that two regions in agreement with all current experimental constraints exist where the 
final state radiation processes are suppressed with respect to the annihilation into photons. Thus, 
there is hope that actually a line could be seen over the continuum background in these regions and 
more information about the dark matter particle could be extracted. We do not discuss the reach 
of proposed future experimental searches in detail but rather focus on general features necessary to 
be able to distinguish lines from the continuum. 

This article is organized as follows. In Section II we discuss the main properties of the simplest 
spinless dark matter model as well as all relevant experimental constraints including the relic density. 
In Section III we discuss in great detail the possible gamma lines in this model and the correlation 
between the gamma rays coming from final state radiation and the annihilation into 77 and Z 7 . 
In Section IV we summarize our main results. In the appendix we list all relevant formulas used in 
this article. 


II. SCALAR SINGLET DARK MATTER 

A. The Model 

In the scalar singlet dark matter model (SDM) the dark matter candidate is a real singlet scalar 
held S which interacts with the Standard Model (SM) particles through the Higgs portal [7]. The 
Lagrangian of this model is very simple and is given by 

/:sDM = /:sM + - XsS^ - \h^hs\ (i) 

where H ~ (1, 2 , 1/2) is the SM Higgs boson and TsM is the usual SM Lagrangian. Once the SM 
Higgs acquires a vacuum expectation value, (H) = vojy/2 where vq = 246.2 GeV, the physical mass 
of the dark matter candidate reads as 

Mg = nig + XpVQ. ( 2 ) 

As is usually done, we assume a discrete Z 2 symmetry to guarantee dark matter stability. Under 
this symmetry S —)• —S such that all odd terms in the scalar potential are forbidden. 

Once the electroweak symmetry is broken, the dark matter candidate S can annihilate into all 
Standard Model particles through the portal coupling Xp. In this model one has only two relevant 
parameters for the dark matter study, the physical dark matter mass Mg and the Higgs portal 
coupling Xp. This is the reason why one can make definite predictions in this model once the relic 
density constraints are used. This model can be considered as a toy model for dark matter, but also 
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is the perfect scenario to understand the possible predictions for different experiments and their 
interplay. 


B. Experimental Constraints 

1. Higgs Decays 

The most conservative, model-independent limit on the Higgs invisible decay branching ratio is 
set by CMS to be BR(/i —)• inv) < 0.58 [26]. However, if we study the predictions for the invisible 
Higgs decay in a particular model, the situation can be rather different. In the scalar singlet dark 
matter model, there is no modification to Higgs physics at the LHC apart from a possibly large 
invisible decay to dark matter if allowed kinematically. Since also the Higgs production cross section 
is unaffected in this model, the invisible width modifies the signal strength of the Higgs decay to a 
P 1 P 2 final state in the following way: 

cj X BR(/i ^ P 1 P 2 ) _ BR{h P 1 P 2 ) 


RPlP2 = 


■nSM 

h 


= 1 — BR(/i —>• inv). (3) 


cjSM X BR{h PiP2)^^ BR(/i ^ PiP2)^^ r™ + r”v 
The combined limit from the final states WW*, ZZ*, 77 , bb, is given by Rtotal = 1-17 ± 

0.17 [27]. This leads to a 95% confidence upper bound on the invisible Higgs branching ratio of 


BR(/i —inv) < 0.16. 


( 4 ) 


This bound is valid for any model which only modifies the Higgs invisible branching ratio. Note that 
a statistically significant deviation of the combined signal strength above one, Rtotal > 1 ) would 
rule out this simple dark matter model up to Mg = Mg/2. The bound obtained here is used when 
we later show the allowed parameter space in the low-mass region, see Figs. 6 and 7 for details. 

2. Relic Density 


In order to compute the relic density of our dark matter candidate S, we use the analytic 
approximation [28] 

1.07 X 109 GeV^ 


= 


( 5 ) 


J{xf)^ Mpi 

where Mpi = 1.22 x 10^® GeV is the Planck scale, 5 * is the total number of effective relativistic 
degrees of freedom at the time of freeze out, and the function J{xf) reads as 

= (o-Urel )(a:) 
fXf ^ 


r 

Ji^f) = / 

J Xf 


-dx. 


( 6 ) 
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Ms [GeV] 

Figure 1. Dark matter relic density VLsh? as a function of the dark matter mass Ms for different values of 
the portal coupling: Ap = 0.1 (solid green), Ap = 0.02 (dashed red), and Ap = 0.001 (dotted blue). The thin 
band where S makes up the full DM relic density today, = 0.1199 ± 0.0022 [29], is marked in light 

blue. 


The freeze-out parameter Xf = Ms/Tf can be computed by solving 


Xf = In 


0.038 g Mpi Ms {avrei)ixf) 

yj g*Xf 


(7) 


where g is the number of degrees of freedom of the dark matter particle. Details on the calculation 
of the cross sections for the different DM annihilation channels and the corresponding analytic 
formulas, including the expressions to perform the thermal average of the cross section times relative 
velocity (crUrei), can be found in Appendix B. 

In Fig. 1 we show the dark matter relic density as a function of the mass Mg, for different 
values of the portal coupling Ap. Depending on the value of the coupling, the correct relic density 
can be achieved in the low-mass regime around the resonance at half of the SM Higgs mass M^ = 
125.7 GeV [27], or off the resonance in the high-mass regime. For some values of Ap there is a 
solution in both regimes. 


In order to understand the different annihilation channels relevant for our study we show in Fig. 2 
the cross sections times velocity crUrei for the different DM annihilation channels in agreement 
with the DM relic density constraints. For every value of the dark matter mass Ms, we use 
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Figure 2. Cross sections times velocity avrei for the relevant dark matter annihilation channels as a function 
of the dark matter mass Ms, setting the coupling Xp such that we have the correct relic density for every 
value of the dark matter mass. 

the corresponding coupling Xp that results in today’s full DM relic density, = 0.1199 ± 

0.0022 [29]. As expected, in the low mass regime the dominant channels are bb and and 

after threshold the annihilation into W~^W~ and ZZ become dominant. Below Ms = 150 GeV, 
we calculated the cross sections from the tabulated partial Higgs widths [30], such that three- and 
four-body decays of the gauge bosons below threshold as well as QCD corrections are included; see 
Appendix B for details. In the high-mass regime, the contributions from the annihilation into the 
SM Higgs h and top quark pairs are significant. 



3. Direct Detection 

To discuss the possible constraints from dark matter direct detection experiments we need to 
know the elastic nucleon-DM cross section. In the scalar singlet DM model, the spin-independent 
nucleon-DM cross section is given by 

_ A 2 /^^ 2 ^^ 
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Ms [GeV] 

Figure 3. Spin-independent nucleon-DM cross section ctsi. The prediction for the cross section shown here 
(black solid curve) is in agreement with the relic density constraints. In blue (dashed) we show the LUX 
bounds [31] and in green (dotted) we show the future reach of XENONIT [32]. In orange (dash-dotted) 
we show the coherent neutrino scattering background [33[. The red part of the curve is excluded by the bb 
limits from Fermi-LAT [34[, see Fig. 5 for more details. 

where misf = (mp-|-mn)/2 = 938.95 MeV is the nucleon mass for direct detection, /jv = 0.30±0.03 
is the matrix element [20], and /r = m]\fMs/{rriM ± Ms) is the reduced nucleon mass. 

In Fig. 3 we show the predictions for the spin-independent nucleon-DM cross section cJsi for the 
typical choice /jy = 0.30 and the corresponding experimental bounds. This model is very simple 
and one can predict clearly the values for the elastic cross section once the relic density constraints 
are imposed. As it is well known, the experimental bounds assume that the dark matter particle 
under study makes up 100% of the DM of the Universe. An important observation is that around 
the Higgs mass resonance direct detection experiments are not able to probe the parameter space 
in the near future as apparent in Fig. 3. However, as will be shown later, indirect searches are 
particularly sensitive to the resonant region and thus highly complementary to direct detection 
experiments. The projected limits by XENONIT [32] tell us that one can test this model for a dark 
matter mass up to a few TeV. 

Note that the current limit from LUX [31] on the scalar singlet DM model strongly depends 
on the particular value that is chosen for /^r. In Fig. 4, we show the predictions for /jv = 
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Ms [GeV] 


Figure 4. Spin-independent nucleon-DM cross section ctsi for different values of /at. In blue (solid) we show 
the current LUX bound [31]. 


0.27, 0.30, 0.33. Depending on /jv, the limit on the dark matter mass Ms varies between 86 GeV 
and 106 GeV. 


4- Missing Energy Searches 

As it is well known, one can hope to observe missing energy signatures at colliders from the 
presence of a dark matter candidate. We have discussed the low mass region where using the 
invisible decay of the Standard Model Higgs one can constrain a small part of the allowed parameter 
space in this model. Unfortunately, in the resonance region the invisible branching ratio of the Higgs 
can be very small and one cannot test this model in the near future. In the heavy mass region one 
can use mono-jets and missing energy searches where one produces the scalar singlet through the 
Standard Model Higgs. Unfortunately, in this case the production cross sections are small and this 
analysis is very challenging. See Ref. [35] for a recent discussion. 
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Figure 5. Velocity-averaged cross sections times velocity for the different relevant indirect detection channels 
using the value of the portal coupling that gives the right relic density. We show the corresponding bounds 
from the Fermi-LAT [34, 36] and H.E.S.S. collaborations [37]. The red parts of the curves are excluded by 
the LUX direct detection limits [31], see Fig. 3 for more details, (a) Annihilation into bb and 77 : the blue 
(solid) line shows the prediction for the annihilation into two b quarks, while the black (dashed) line shows 
the prediction for the annihilation into 77 . The dip and subsequent increase in the DM annihilation cross 
section to 77 is a feature of the loop functions involved in h* —>■ 77 , see Eqs. (B10)-(B13) in Appendix B for 
more details, (b) Annihilation into bb and Z 7 : the blue (solid) line shows the prediction for the annihilation 
into two b quarks, while the black (dashed) line shows the prediction for the annihilation into Zj. 


5. Indirect Detection 

In Fig. 5 (a) and (b) we show the predictions for the dark matter annihilation into bb as well as 
77 and Z'j, respectively, together with the bounds from the Fermi-LAT [34, 36] and H.E.S.S. [37] 
experiments. Continuum searches in bb and line searches constrain the same mass region, with the 
continuum searches giving the slightly more restrictive bounds. These bounds are very important 
because one can rule out part of the parameter space close to the resonance region. This is the 
only way to exclude this region because the contribution to the invisible decay of the Higgs is very 
small. Unfortunately, in the heavy mass region the current experimental bounds cannot exclude 
any of the parameter space. However, in the near future these experiments could test this simple 
model if the dark matter mass is close to 100 GeV. 
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Ms [GeV] 

Figure 6. Allowed parameter space in the Mg-Xp plane in agreement with the relic density constraints, 
direct and indirect detection, and invisible Higgs decays. In gray we show the region of the parameter 
space where one overdoses the Universe; the black line corresponds to today’s full relic density, = 

0.1199 ± 0.0022 [29]. In green we show the bounds from the invisible decay of the SM Higgs, using the 
CMS bound BR(/i —>■ SS) < 58% [26], as well as the calculated limit from Eq. (4). The red part of the 
relic density curve is excluded by the LUX direct detection experiment [31], while the blue part of the curve 
shows the projected reach of the XENONIT experiment [32]. The orange part of the curve is excluded by 
the hh limits from Eermi-LAT [34]. 


6. Summary 

In Fig. 6 we show the allowed parameter space in agreement with the relic density constraints, 
direct and indirect detection, as well as invisible Higgs decays. As one can appreciate, there are two 
main regions allowed by all experiments. In the low mass region, 53 GeV < Ms < 62.8 GeV, the 
dark matter annihilates through the Higgs resonance, while in the heavy mass region, Ms > 92 GeV, 
all the gauge boson channels are open and dominate. In Fig. 6 we also show the experimental 
bounds on the invisible decay of the SM Higgs [26] and the projected direct detection bounds from 
the XENONIT experiment [32]. The gray region is ruled out by the relic density constraints because 
in this region one overdoses the Universe having too much dark matter relic density. Notice that 
even this simple model for dark matter is not very constrained by the experiments. 
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Figure 7. Allowed parameter space in the Ms-Xp plane in agreement with the relic density constraints, direct 
and indirect detection, and invisible Higgs decays in the low mass regime. In this region the constraints 
from the invisible decay of the Higgs shown in green are very important. Color coding is the same as in 
Fig. 6. 

Low Mass Regime 

In the low mass region the allowed dark matter mass is 53GeV < Ms < 62.8 GeV. In this 
region close to the Higgs resonance the dark matter can annihilate into Standard Model fermions 
or into two fermions and a gauge boson. In Fig. 7 we show a detailed analysis of this region to 
understand which part of the parameter space is ruled out by experiments. Notice that the main 
annihilation channel is SS —)• bb. In this model one can set bounds only using the constraints on 
the nucleon-DM cross section. The scattering between electrons and DM is highly suppressed by 
the small Yukawa coupling. 

From the results presented in Fig. 7 one can see that the resonance region cannot be excluded 
or tested in the near future by direct detection experiments. This is a pessimistic result, but 
fortunately this region can be tested at gamma-ray telescopes as we will discuss in the next section. 

Heavy Mass Regime 

When the dark matter is heavy it can annihilate into all Standard Model particles. For 
Ms ^ 70 GeV, where the dominant annihilation channels are into gauge bosons, the current di- 
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rect detection bound from LUX [31] rules out only a very small range in the parameter space, 
Ms ^ 92 GeV, see Fig. 6 . It means that this model is not very constrained by direct detection in 
the heavy region. These results are crucial to understand the testability of this simple model at 
gamma-ray telescopes. 


III. DARK MATTER ANNIHILATION INTO GAMMA-RAY LINES 

The detection of a monochromatic gamma line coming from dark matter annihilation would be 
a very strong hint towards the particle interpretation of the dark matter in the Universe. It is 
highly unlikely that an astrophysical compact source would generate very energetic monochromatic 
photons in different regions of the galaxy. The question whether a gamma line signal is generic and 
visible in a dark matter model is very subtle. Whenever tree-level annihilations into SM particles 
are present, one expects that the final state radiation off the charged SM particles can make the 
gamma line undetectable experimentally. To investigate the line visibility in the scalar singlet dark 
matter model we compute the cross sections for the final state radiation processes and investigate 
the gamma flux spectra in the relevant parts of the parameter space. 

A. Final State Radiation 

There are three relevant regions which define the properties of the gamma spectrum coming 
from dark matter annihilation. Let us define = E^/Ms, where is the energy of the photon 
and Ms is the dark matter mass. When is very small one has the photons coming mainly 
from hadronization, i.e., the dark matter annihilates into quarks and from the cascade one has 
the photons with a continuous spectrum. When Xry is close to one, one finds that the final state 
radiation processes contribute more because they can provide hard photons. Finally, when x^ = 1 
one has the gamma line with energy equal to the DM mass. Therefore, one must understand the 
hnal state radiation processes to investigate the visibility of the gamma lines. 

The relevant final state radiation process for our study is SS —>• XX^, with the kinematic 
endpoint of the continuous 7 spectrum at 

= Ms (1 - w) . (9) 

in the non-relativistic limit s = 4M|. In the low-mass regime the dominant process is SS —)■ //y 
with the strongest contribution from the bottom quark, while in the high-mass regime SS —>• 
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W~^W 7 becomes dominant. The differential cross section times velocity of those processes is 
given by 


1 


dcrVrel _ _ 

dE^dEi 327r3s 


|A 4 FSRr 


where the integration limits for the integration over Ei for a fixed E^ are given by 


= Ms- 


^max = Ms - 


+ 


^E^{E^ - Ms)Ms [Ml + {E^ - Ms)Ms] 

2Ms{Ms - E^) 

^E^{E^ - Ms)Ms [Ml + {E^ - Ms)Ms] 

2Ms{Ms - E^) 


( 10 ) 


( 11 ) 

( 12 ) 


in the limit s = 4Mj. See Appendix C for the amplitudes of the two relevant processes for final 
state radiation, SS —)■ ffj and SS —)• Notice that in the low-mass region the FSR is 

suppressed by small Yukawa couplings. Therefore, this is the region where generically one can have 
a visible gamma line. 


B. Gamma Flux 


The differential photon flux is given by 

d^-y Ury d{avrei) j / 

^ ^ dE^ ^ ’ 

where the factor Jann contains the astrophysical assumptions about the DM distribution in the 
galaxy and thus all the astrophysical uncertainties. Here is the number of photons per annihi¬ 
lation, and dN^/dE^ is the differential energy spectrum of the photons coming from dark matter 
annihilation. In all numerical calculations, we will use the J-factor from the R3 region-of-interest, 
given by the Fermi-LAT collaboration to be Jann = 13.9 X 10^^ GeV^cm“^ [38]. The R3 region 
is a circular region of radius 3° centered on the galactic center [38]. The differential flux of the 
line is extremely narrow, however to make connection with the experiment it will be folded with a 
Gaussian function modeling the detector resolution. 



• SS —)• 77 : for the annihilation into two photons the flux is given by 


dE^ 


e/w;, E^) , 


with 


exp 

V^Eo{^/w) 


(14) 


G{E^,i/w,Eo) 


(15) 
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Table 1. Benchmark scenarios for the study of the 7 spectrum. 


Scenario Ms [GeV] 

Xp 

Energy of the Z 7 line [GeV] 

1 

62.5 

9.06 X lO-'^ 

29.2 

2 

150 

2.08 X 10-2 

136 

3 

316 

4.17 X 10-2 

309 

4 

500 

6.87 X 10-2 

496 


The parameter ^ is a measure of the detector energy resolution which varies between 0.01 
and 0.1 in the relevant energy range. The factor w = 2-^/2 log 2 ~ 2.35 determines the full 
width at half maximum as aow = ^Eq, therefore we have (Tq = Eq^/w in the usual Gaussian 
function. For the annihilation to 77 , the energy of the gamma line is at the dark matter 
mass, 


E'y = Ms- 


(16) 


SS —)• Xj: for the annihilation into an unstable final state particle along with a photon, the 
flux is given by 


Jann {o'Vrel)X'f 


dE^ 


i 


00 

dEo- 


AMsMxTx 


^G{E^,i/w,Eo). (17) 


SvrMl Jq tt {4.MI- AMsEq-+ 

Here Tx is the decay width of the unstable particle in the final state and Mx is its mass. 
See Appendix D for a derivation of the differential energy spectrum used in Eq. (17). The 
gamma line energy is given by 

Ml 




(18) 


Using these expressions for the differential flux we now study the predictions for the gamma lines 
in this model in the benchmark scenarios defined in Table I. 

In Fig. 8 (a) we show the gamma spectrum for the scenario 1 with Ms = 62.5 GeV. In this 
case one has the resonant dark matter annihilation through the SM Higgs. As one can see in this 
scenario it is possible to identify the gamma line from DM annihilation into 77 , while the line from 
Zj is not visible in the plot since it is at x-f = 0.47 and will be swamped in the FSR background. 
The main contribution to final state radiation in this case is coming from the annihilation into 667 
but it is suppressed by the small bottom Yukawa coupling. Therefore, in this case one has a large 
difference between the final state radiation and the gamma line. 
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Benchmark Scenario 1: Ms = 62.5 GeV, \p = 9.06 x 10 ^ 



Benchmark Scenario 3: A/g = 316 GeV, Ap = 4.17 x 10 ^ 



x~i = E^/Ms 

(c) 


Benchmark Scenario 2: A/g = 150 GeV, Ap = 2.08 x 10 ^ 



Xj = E^/Ms 

(b) 

Benchmark Scenario 4: A/g = 500 GeV, Ap = 6.87 x 10“^ 



x~^ = E-i/Ms 

(d) 


Figure 8. Spectra for the benchmark scenarios from Tab. I. The different curves correspond to an energy 
resolution of 1% (solid green), 5% (dashed red) and 10% (dotted blue), respectively, (a) Benchmark scenario 
1: the main contribution to final state radiation is the annihilation to bb-j. (b) Benchmark scenario 2: the 
main contribution to final state radiation is the annihilation to WW^. (c) Benchmark scenario 3: the main 
contribution to final state radiation is the annihilation to WW^. (d) Benchmark scenario 4: the main 
contribution to final state radiation is the annihilation to WW^. 


In Fig. 8 (b) the predictions for the gamma spectrum is shown for the second scenario where 
Ms = 150 GeV. This scenario is ideal because one can see the two possible lines in this model, the 
77 and Z 7 lines, if one has a good energy resolution. In this case the main contribution to hnal 
state radiation comes from the DM annihilation to WW^. However, as one can appreciate, there 
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Ms [GeV] 


Ms [GeV] 


(a) 


(b) 


Figure 9. Line visibility for the three different experimental energy resolutions used before: 1% (solid green), 
5% (dashed red) and 10% (dotted blue). To estimate the visibility of the line, we mark the ratio 1 in black 
(dash-dotted) and the ratio 10 in orange (dash-dotted), (a) Visibility of the 77 line: we show the ratio of 
the differential gamma flux at the line energy to the FSR flux at 90 % of the dark matter mass for 10 % 
energy resolution, at 95 % for 5 % energy resolntion, and at 99 % for 1 % energy resolution, (b) Visibility 
of the Zj line: we show the ratio of the differential gamma flux at the line energy to the FSR flux at 90 % 
of the line energy for 10 % energy resolution, at 95 % for 5 % energy resolution, and at 99 % for 1 % energy 
resolntion. 


is a large difference between FSR and the gamma lines because the endpoint of the FSR is far from 
the DM mass. 

The case when the DM mass is 316 GeV is shown in Fig. 8 (c). There is a large difference 
between the rate for the Z 7 and 77 lines. Unfortunately, in this case one could see the lines only 
with a perfect energy resolution. The cross section for the hnal state radiation processes is large in 
this case making the observation of gamma lines very challenging. Finally, we present in Fig. 8 (d) 
the energy spectrum for the case when Ms = 500 GeV. In this case one cannot distinguish the 
gamma lines since the difference between the final state radiation and the gamma line is very small. 

In order to have a more generic discussion about the visibility of the gamma lines in Fig. 9 we 
show the ratios between the 77 and Zj fluxes and the gamma flux from final state radiation. We 
display these ratios for Ms > 100 GeV. We show the curves for 1% energy resolution (green solid), 
5% energy resolution (red dashed), and 10% energy resolution (blue dotted). To be conservative 
we can say that the lines are visible if the ratio between the fluxes in Fig. 9 is larger than a factor 
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10. This means that for realistic experiments with an energy resolution of 5 %, the gamma lines 
can be visible when the dark matter mass is smaller than 300 GeV. 

As one can appreciate from the above discussion one could observe the gamma lines from dark 
matter annihilation in this model only when the dark matter mass is small, i.e., in the low mass 
region where Mg = (53 — 62.8) GeV or in the intermediate region, Mg = (92 — 300) GeV, where the 
two gamma lines could be distinguished from the continuum and among each other. These results 
are crucial to understand the testability of this model at gamma-ray telescopes. 

In Fig. 5 one can see that the Fermi-LAT experiment could test this model in the low mass region 
where the gamma gamma cross sections are large. The proposed GAMMA-400 experiment [39] will 
have a better energy resolution and therefore will explore the same region with a better chance to 
see the gamma lines. For the high-mass region, there is currently no planned experiment that will 
have the required sensitivity in cross section to reach the values predicted in this model. 

Now, let us make a short summary of all constraints to understand how this model can be tested 
in the near future. In Fig. 3 we have seen that the low mass region cannot fully be tested at direct 
detection experiments. However, the region where the dark matter mass is in the range Mg = 
(92 — 300) GeV could be tested at the XENONIT experiment. Notice that the spin-independent 
cross section for Mg = 100 GeV is crsi = 8.6 x 10““^® cm^, while for Mg = 300 GeV one finds 
cjsi = 6 X 10“^® cm^. At the same time in this region, one could see the gamma lines coming from 
dark matter annihilation as we have shown in Fig. 9. The annihilation cross sections are ((TUrel )77 = 
9.1 X 10“^^ cm^s“^ when Mg = 100 GeV, and (crurei)z 7 = 5.2 x 10“^^cm^s“^ for Mg = 300 GeV. 
Therefore, one could say that only the intermediate region with Mg = (92 — 300) GeV can be tested 
at both direct and indirect experiments. Of course, due to the smallness of the gamma line cross 
sections, it will be challenging to reach the required experimental sensitivity at indirect detection 
experiments to test this model at Mg = 300 GeV. 

IV. SUMMARY 

In this article we have investigated in great detail the predictions in the simplest dark matter 
model where the dark matter candidate is a real scalar field. This model can be considered as a 
toy model for dark matter but it offers the possibility to connect the predictions for all relevant 
dark matter experiments. This model has only two free parameters which are constrained by the 
relic density. Once one uses the relic density constraints we can make predictions for the elastic 
nucleon-DM scattering and the cross sections for DM annihilation into gamma rays. We have 
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revisited the model and updated all constraints for the full parameter space of the model. 

We have computed for the first time all contributions to final state radiation analytically in this 
model. In the low mass region, Ms = (53 — 62.8) GeV, the main contribution to final state radiation 
comes from the annihilation into two 6 -quarks and a photon, while in the heavy mass region the 
annihilation into two W gauge bosons and a photon is the dominant contribution. These results 
are very important to understand the visibility of the gamma lines. 

We have shown the predictions for the two possible gamma lines in this model coming from 
the DM annihilation into two photons and into Z'j. In Section III we have shown the numerical 
predictions for the gamma spectrum in different scenarios. We have shown that in the low mass 
region one could see the 77 line because the annihilation into 667 is suppressed. In the intermediate 
region where the dark matter mass is Ms = (92 — 300) GeV there is a large difference between the 
continuum and the gamma line and one could see the two gamma lines if the energy resolution is 
good. Unfortunately, in the heavy mass region the dark matter annihilation into IUVU 7 is large 
and it is very challenging to observe the gamma lines. We have shown that only the region where 
Ms = (92 — 300) GeV can be tested at direct and indirect detection experiments. These results 
can motivate a dedicated study in this region of the parameter space if one sticks to the simplest 
model for dark matter. 
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Appendix A: Invisible Higgs Decay 

In the scalar singlet DM model the Standard Model Higgs can decay into dark matter in the 
low mass region. The invisible decay width of the Higgs in this case is given by 

r(h ^ SS) = - 4M|)1/2. (A1) 

To calculate the invisible branching fraction 

T{h^ SS) 

PSM + r(/i ^ SS) ’ 


BR{h SS) 


(A2) 
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we use the SM Higgs width F™ = 4.17MeV for = 125.7GeV [30]. 


Appendix B: Dark Matter Annihilation Cross Sections 


• Annihilation cross section times velocity from the Higgs width: to obtain an accurate value 
for the total DM annihilation cross section to SM final states, one needs to take into account 
QCD corrections for quarks in the final state, as well as three- and four-body final states from 
virtual gauge boson decays below threshold. This is most easily done by using the tabulated 
Higgs width^ as a function of the invariant mass [30] and rewrite the cross section 

times velocity as 

\Dh{s)f rh{V~s), (Bl) 

Vs 

with 


(s - M2)2 + MlTliMh) ■ 

This factorization is possible for all final states except the SM Higgs, such that for Ms > 
the contribution SS —)• hh has to be added. For Ms < Mh/2, the width in Dh{s) has to take 
into account the invisible decay h —)• SS. The thermally averaged annihilation cross section 
times velocity (crurei) is a function of x = Ms/T, and - when using the expression for av^ei 
in Eq. (Bl) - can be computed as 

(iTUrel)(x) = - e-:t—— [ ds\/s — 4Mn S Ki ( CUrel, (B3) 

where Ki and K 2 are modified Bessel functions of the second kind. 


For Ms > 150 GeV we use the tree-level expressions calculated below, since then loop cor¬ 
rections overestimate the tabulated width. See the discussion in Ref. [20] for more details. 
From the cross section cr, the thermal average can be computed via 

(o-Urei)(x) = [ ds {s - AMI) ^fs Ki (a. (B4) 

8M^K^{x) \Ms J 

• Annihilation into Standard Model fermions: 


a{SS ^ //) 


XIM^nI{s - 4M2)3/2 


27:8x1 s — 4M| 


{s-Mlf + TlMl 


(B5) 


where nI is the color factor of the fermion /. 


^ Or the partial widths for each final state when we are interested in the individual cross sections times velocity. 
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Annihilation into two W gauge bosons: 

A 2 ^s - 4M^ (s 2 4 ^^^ + i 2 M^) 


^ lyw^) = 


4^" - 4M| [(.-M^^ + r^M^ 

Annihilation into two Z gauge bosons: 


a(SS ZZ) = 


Annihilation into two Higgs bosons: 


Ap ^ 

/s - 4M| 

(s2 - 4M|s + 12M|) 

87rs ^ 

/»- “VJ 

/»- vj)' + rjM| 




16AH 


•p'^O 


y/s - 4M| [. ^h-4M|M| + M|s 

32Apt>g (s^ - 2ApSt>^ - 4M^ + 2ApM|i;^) 


2 M 2 - s 


[2Mt - 3Mls + 

In the limit s —)• 4:Mg, the cross section times velocity is given by 

.„„.(SS ^ HH) = A;W-4Mj + 2Vg(4Mj-Ma]^ ^^ 

47rM| (M^ - 6M2 m| + 8M|) V 


Annihilation into 77 : 

a{SS —^ 77 ) = 

where the width is given by [40] 


02.3/2 

r (s) = ---^ 

^ 2567r3?;2 


4Ap Vq r^^(s) 


s-4M| 




^NiQ}A>l/,{rf) + A^,{rw) 
f 


with the form factors 


Al,,{T)=2[T + {T-l)f{T)]T-\ 

AJ'(t) ^ - [2t^ + 3t + 3(2r - l)/(r)] T~^, 


and the function 


/(t) = ^ ^ _ 

The parameters r* for fermions and the W are given by 


arcsin^yT 


12 


r < 1 


r > 1 


5 S 

Tf = 777 , and t\y = 


AMj 


4M^' 


(B 6 ) 


(B7) 


(B 8 ) 


(B9) 


(BIO) 


(Bll) 

(B12) 

(B13) 


(B14) 


(B15) 










































22 


Annihilation into Z 7 : 


a{SS Zj) = 


4Ap Vo rz^(s) 


s - 4A/| 




(B16) 


The width is given by [40] 






with 


'^N;[^^Ay^{Tf,Xf) + A’^irwAw) 
^ cw ' 


(B17) 


Vf — 21 ^ — AQfs'^, 


(B18) 


and the parameters^ 


4M2 ^ ^ 4M2 

Ti = - and Xi = — 

s Mi- 


(B19) 


The form factors are 


^ 1/2 = [h h (r, A)] , 

"1^ ~ c^) -^ 2 ( 7 ', A) + 





/i(r,A) , 


with the functions 

= 2 (V^ + ^ 

^2(^, = ~ 2{t-X) ~ ’ 

with /(r) from Eq. (B14) and 


r^A 

(r - A )2 


[ai-r -5(A ^)] , 


(B20) 

(B21) 


(B22) 

(B23) 


9{r) 


7 

1 

1 

1 arcsinyT" 

Vl—T ^ 

2 



r > 1 
r < 1 


(B24) 


Appendix C: Final State Radiation 


In this appendix, we use 



and 


= 


Ws' 


(Cl) 


Notice that the definition of the Ti is the inverse of the definition used in the /i —>■ 77 case. 


2 
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SS —)• 7 / 7 : in the limit s —)• 4Me the the matrix element used in Eq. (10) is given by 


sxyqjMjNi 


|a^fsr(55 ^ /77)|' =- 

Mlixi - l) 2 (xi +x^- 1)2 ^(M 2 - 4M2)" + r 2 M 2 

X |m| (xi — 1)^ + MjMg \^x\ + 2xf (xj — 3) — 2xi (x-y — 3) Xy + Xy (5xy — 8) + 2 ] 

— 2Mg (xi — 1) \{xy — 2)xy + 2] (xi + — 1) |. (C2) 

Integrated over xi, we obtain for the differential cross section 


daVre\{SS 7 / 7 ) 


xyqjMjNi 


dXr, 


dvr^M^x^ 


(M2-4M2)Vr2M2 


X I (M| + 2M|) y {xy - 1) [m2 + Ml{xy - 1) 

( 

+Ms \2M‘j{xy — 1) + M| (2 + Xy{xy — 2))] In 


1 + ■»/1 + 


\ 




1 -A 1 + 




(^7 1)-^S / , 


7 (C3) 


S'S —^ IEVE 7 : in the limit s —^ 4Mj the matrix element used in Eq. (10) is given by 


|A^FSR(55^1ElE7)r = 


4A2e2 


(xi - if (xi + x.^ - if (4M2 - m 2)" + M2r2 


X I — 3M^Xy + 16M| (xi — 1) (x.^ — 1) (x-y + xi — 1) + 4M|M^ [4 + dx" (1 + (x-y — l)x-y) 
+ 4xi (Xy, — 2) (1 + (Xy — l)Xy) + Xy ( —8 + Xy(7 + 2(Xy — 2)Xy))] 

+ 4M|M^ [ — 3 + 3Xy (Xy — 1) + 3xi (Xy — 2) (Xy — 1) — 2 (Xy — 3) Xy] |. (C4) 


Integrated over xi, we obtain for the differential cross section 


daVreijSS WW-i) 
dX-y 


e^Xl 


STT^MgXy 


(ml - mIY + Afjrj 


X - 1) [AfJ, + MUx., - 1)] [mlMl, - 3Mi. + 4AfJ (2x1 - 4] 

(l- ~~ 

— (4MJ — dMjM-^ + 3M^) [M^ + 2Mj(xy — 1)] In - ^ . (C5) 


I — \ I + 




(^7 ^)^'s ) , 
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Appendix D: Gamma-Ray Spectrum from Dark Matter Annihilation 


We derive the expression for the differential cross section of dark matter annihilation to a photon 
and an unstable particle X. Consider the annihilation SS —)• X -|- 7 —)• ff'j. We begin by 
decomposing the three-body phase space of the final states into two two-body phase space parts 
and an integral over the mediator as 


dU{y/s]pi,p2,p-y) = 


d?pi cfp2 


d?p^ 


{2'Kf2Ei {2'Kf2E2 {2'Kf2E. 


■{2-K)'^5{ki + k2-pi-p2- P'f) 


d^p d^p^ 


-{27r)^6{ki +k2 - p-p^)- 


d^pi d^p2 


-{2'nf5{p-pi -P 2 ), 


(Dl) 


(27r)4 (27r)32£;^' ' ' (27r)32£;i (27r)32^2 

where pi and p 2 are the momenta of / and /, p^ is the photon momentum, p the off-shell momentum 
of X and the invariant mass of pi and p 2 combined, fci and ^2 are the momenta of the incoming 

one finds that 


particles. Using the relations = j^{Pq — |pp) = —2|p| and EpdEp = 

d^Py 


dM‘^ = —2EpdEp and 


dU{y/s;pi,p2,Py) = - 


d\p\ d\p\ '”0 

dM"^ d?p 


2ti {2tt)^2Ep {27r)^2E- 


-{27r) 6{ki + k2-p-py) 


X 


d^pi d^p2 


-{2tt)'^6{p - Pi -P 2 ) 


{2'Kf2Ei {2'Kf2E2 
= — ^dU{^/s■,p,py)dU[p;pl,p2). 

With this decomposition, the cross section can be written as 

|A4P \M\‘^ dM'^ 

da = — — d^{^/s,pl,p 2 ,Py) = -r— dU{^/s■,p,py) dU{p]pi,p 2 ), 

<Pss 5*55 27r 

where $55 is the dark matter flux. Now, we write the differential amplitude of the process assuming 
a narrow width of the intermediate particle. Then, one finds 


(D2) 


(D3) 


da 

dM^ 


1 


1 


— m\ + iVxmx 


[ dIl{^/s■p,Py) f \Mx^qq\‘^dU{p;pi,p2) 

Jp,p^ ^SS •Jpi,P2 


a{SS —)• X 7 ) Txmx 


IT 


1 


+ iTxmx 


(D4) 


Using the fact that = 4M| — AMsEy, which gives a Jacobian factor of = —AMs, we obtain 
the final result for the differential cross section 

da a{SS —)• X 7 ) AMsmxTx 


dE^ 


vr (4M| - 4 M 5 £;.^ - m^)2 + 


^aiSS^Xj)^. (D5) 
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